Lifetime and Emission Characteristics of Electronic Excitations in 2D Optical Lattices 
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Collective electronic excitations "excitons" in planar optical lattices exhibit strong modifications 
of the radiative damping rate and directional emission pattern as compared to a single excited 
atom. Excitons for long wave numbers and polarizations orthogonal to the lattice plane exhibit 
superradiance with a very short life time and a tightly confined emission direction. For shorter 
wavelength and in plane polarization they can posses a long life time, which beyond a critical wave 
number tends to infinity. Those excitons thus become metastable and decoupled from the free 
radiation field. They can store a single energy quantum for a long time and transfer the excitation 
over long distances. In general the spatial, polarization and frequency dependence of the emission 
pattern can provide us with optical and electronic properties of optical lattices. 

PACS numbers: 37.10.Jk, 42.50.Pq 



Optical lattices with ultracold atoms constitute an es- 
tablished test system for a wide range of condensed mat- 
ter models and phenomena i ij, l2| . Recent theoretical and 
experimental progress highlights the importance of such 
a system both for fundamental physics and applications. 
Ultracold atoms loaded to optical lattices are well de- 
scribed by the Bose-Hubbard model Q that predicts the 
quantum phase transition from the superfluid into the 
Mott insulator phase [4] . The big similarity between op- 
tical lattice ultracold atoms in the Mott insulator phase 
and solid crystals encouraged us to study different solid 
state effects for such a set-up. Many advantages of opti- 
cal lattices over solid crystals are a result of the precise 
controllability of the system properties, e.g. the number 
of atoms per site, the geometry and symmetry. These 
advantages opened the door for a deep understanding of 
different solid state effects and to answer open questions 
in the field, beside the emergence of new phenomena 

At very low temperature in optical lattices the atoms 
are localized in the lowest Bloch band forming a Mott 
insulator phase. The atoms retain their identity except 
from light shifts due to the external laser fields that form 
the optical lattice [1, 0] . This fact makes optical lattices 
analogous to molecular crystals. Optical properties of 
molecular crystals are mainly dominated by the forma- 
tion of collective electronic excitations, in which a single 
molecule electronic excitation is delocalized in the crys- 
tal due to electrostatic interactions and represented by a 
wave that propagates in the material with a fixed quasi- 
momentum and effective mass. Such quasi-particles are 
termed as Frenkel excitons @ . Excitons in optical lattice 
ultracold atoms studied extensively by us [y, 0] ■ 

In the present letter we investigate the radiative damp- 
ing rate and emission pattern of electronic excitations in 
planar optical lattices. Optical excitations in an ultra- 
cold gas are short lived and lead to strong heating and 
a fast decay of coherence via emission recoil and reab- 
sorption, hence most experiments are performed very far 
from resonance. The formation of excitons in 2D optical 
lattices can significantly alter the radiative damping rate 
into free space and its impact on the lattice atoms. We 



show that excitons will decay collectively in superradi- 
ant or subradiant fashion, and can be even metastable. 
As the decay and absorption occur collectively, the re- 
coil and energy difference are spread over large spatial 
regions, depositing only a small amount of energy and 
momentum per particle. A 2D Mott insulator thus can 
only be excited under special conditions and will remain 
robust under near resonant excitations. 

Historically the modified radiative damping rate of ex- 
citons at low dimensional molecular crystals was pre- 
dicted by Agranovich Q , where superradiant decay fore- 
seen for small wave vector excitons. While in those cases 
experimental confirmation are difficult and superimposed 
by many other effects, the situation in cold atom optical 
lattices is much more promising. The special case of ID 
optical lattice was treated by us in [9|]. Experimentally 
the fabrication of 2D optical lattices in the Mott insula- 
tor phase was achieved experimentally by several groups 
. Recently the fluorescence imaging of 2D optical lat- 
tice in the Mott insulator phase is obtained down to a 
resolution of a single site [llj . The emission pattern for 
long wavelength excitons we derive here can serve as an 
observation tool for the kind of excitations exist in the 
system, and to predict the formation of excitons, and 
allow us to extract different properties of the system. 

As our model system we consider a two-dimensional 
planar optical lattice in the Mott insulator phase with 
one atom per site. The optical lattice has a square sym- 
metry of lattice constant a, see figure (1). The atoms 
are taken to be two-level systems of transition energy 
Ea- An electronic excitation can transfer among atoms 
at different sites due to electrostatic interactions with 
the coupling parameter Jnm, where n is the site po- 
sition in the lattice plane. Here n = {nx,ny) with 
n, = a(0,±l,--- ,±N,/2), and N = x Ny is the 
number of lattice sites. The electronic excitations are 
described by the Hamiltonian 



Hex — ^ Ea B^^Bn 



Jnm B^Bin 
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Bl and B„ are the creation and annihilation opera- 
tors of an electronic excitation at site n, respectively. 
The operators are assumed to obey Bose commutation 
relations, which is a good approximation at low den- 
sity of electronic excitations in the lattice. The tran- 
sition dipole operator of an electronic excitation is fi = 

An electronic excitation is delocalized in the lattice 
due to, e.g., resonance dipole-dipole interactions, and 
in exploiting the lattice symmetry it represented by 
a wave that propagates in the lattice with in-plane 
wave vector k. Such collective electronic excitations are 
called excitons. In using the Fourier transform i?„ = 
e*'''"-Bk, we get the diagonal exciton Hamilto- 

nian Hex = J^k^e^C^) bI^Bi^, with the energy disper- 
sion £ea;(k) = Ea + Er •^(R-)e''''^- The coupling pa- 
rameter is a function of the distance between the sites, 
namely R = n — m. For the case of nearest neigh- 
bor interactions with the coupling parameter J, we get 
iJex(k) = Ea — 2 J [cos(fc2;a) -|- cos(fcyo)]. The wave num- 
bers take the values h = ^ (0, ±1, • • • , ±A^,/2). 

The optical lattice is located in free space. The ra- 
diation field of free space is represented by the Hamil- 

tonian Hrad = J2q\^phiq) OqA'^qA. where a]^^ and a^x 
are the creation and annihilation operators of a photon 
with wave vector q and polarization mode A, respec- 
tively. The photon dispersion has the linear dispersion 
Ephil) = hwr{q) = hcq- The electric field operator is 



E(r) 



qA 



flqA eqAe 



zq-r 



4a ^xe-""-] 
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where GqA is the photon polarization unit vector, and V 
is the normalization volume. 

The matter-field coupling is given by the electric dipole 
interaction Hj = fi ■ E, and in the rotating wave ap- 
proximation for linear polarizations we get 



Hi 



qA,n 



hcq 
2^ 



X {a^xBie'^'^-al^^B^e-^^--}. (3) 



We define q = (q', g^), in using the above transformation 
and the relation e^^'^ -k).n _ §^,^^ obtain 



Hi = Y,^fi9^x{b.ixBl-bl^Bu} 

qA 



(4) 



with the coupling parameter ^i^qA 
where now q = (k, (j^), and q'^ = 



hcqN / \ 
2^ ' ^qAj) 

k'^ + q'^, see figure 
(1). The coupling is between photons with wave vector 
component parallel to the lattice equal to the exciton 
in-plane wave vector. 

Here wc calculate the radiative damping rate into 
free space of a single exciton with wave vector k in a 




FIG. 1: The optical lattice plane, with lattice constant a. The 
directions of the in-plane wave vector k and the transition 
dipole n are seen. 



planar optical lattice. We use the Fermi Golden rule 

Tfe = ^j:,^x\{.f\Hi\t)\^S{EeA^)-Eph{q)). The ini- 
tial state is of a single exciton and empty radiation field, 
K) = I lex(k), Op/(), and the final state is of a single photon 
and ground state optical lattice, |/) = lOe^, lp/i(q))- The 

matrix element is given by {f\Hi\i) = l^lv (/^ ' 
The summation over the photon polarization is obtained 
by the relation |/x • eqA|^ = \nf — ^^]^^ ■ The sum 
over is converted into the integral — )• ^ Jq°° dqz, 
where V = SL with S = No?. We define /x ~ (/X|| , /Xx), 
see figure (1). The angle between the in- plane transi- 
tion dipole and the exciton wave number is (j), namely 

^k • = /Z||fccos(^. Furthermore, we define 6 the an- 
gle between the transition dipole and the z direction, 
hence /^y = iisiwO and fiz = HCOsO. Straightforward 
integration yields 



re.(k) 



2 J. El{k) 
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ElA^) 
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EU^) 



X |sin^6»^l- 
- 2 sin ^ cos 9 cos <^-g|I ^£;|,(k) - i?2(A=)| , (5) 

where Eo{k) = tick. The most important result here is 

that for k > kc we have Fk = 0, where hckc = EexO<-c), as 
beyond the singularity point the damping rate becomes 
imaginary. The results need to be compared with the 
radiative damping rate of a single atom in free space, 

which is Tat = 3^3^- 

Here we use typical numbers for optical lattice ultra- 
cold atoms. The lattice constant is a = 1000 A, and the 
transition dipole is = 1 eA. The exciton energy is taken 
to be i?ea;(k) = 1 eV , where for simplicity we neglect the 
k dependence, which is negligible relative to the Eo{k) 
one. We plot the scaled damping rate, Fex(k)/Fat, as a 
function of different directions of the transition dipole, 9, 
different directions of the exciton, (j), and different exciton 
wave vectors, -Eo(^) = tick. 

In figure (2) we plot T^x/Tat vs. Eo{k) for 6 = 0, where 
the transition dipole has only component normal to the 
lattice plane. Here TgxO^) is (f> independent. The damp- 
ing rate starts from zero, where the long wave length 
excitons are metastable, and increases with k to become 
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supcrradiant with damping rate larger that Tat- Close 
to hckc = 1 eV the damping rate diverges, and beyond 
hckc = 1 eV jnmps back to zero, where all excitons be- 
tween kc and the Brillouin zone boundary at n/a are 
metastable, where hcw/a ~ 6.2 eV. 
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FIG. 2: The scaled damping rate F^x/Tat vs. Eo{k) = hck, 
for ^ = 0, and which is independent. 

In figures (3) we plot Tex /Fat vs. Ef){k) for 9 = 7r/4, 
where the transition dipole has equal in-plane and normal 
components. In figure (3a) we take the in-plane wave 
vector k to be parallel to the in-plane transition dipole, 
that is (f) = Q. Very interesting behavior appears here for 
the damping rate. For small in-plane wave vectors the 
excitons are superradiant with damping rate larger than 
Tat, and decreases in increasing k till it becomes zero 
at Eo{k) = i?ej;(k)/\/2 where excitons around this point 
are metastable. For larger k the damping rate start to 
increase and diverges at Eo{kc), and then jump back to 
zero beyond this point. In figure (3b) we take the in-plane 
wave vector k to be normal to the in-plane transition 
dipole, that is (p = ■k/2. Now the small wave vector 
excitons are superradiant and the damping rate increase 
with k to diverge at kc, and then jump to zero. 
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FIG. 3: The scaled damping rate Tex/Tat vs. 
for 6* = 7r/4 with (a) = 0, and (b) </- = 7r/2. 
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Eo{k) = hck, 



In figure (4) we plot Te^/^at vs. Eo{k) for 9 = tt/2, 
where the transition dipole has only in-plane component. 
We take the in-plane wave vector k to be parallel to the 
in-plane transition dipole, that is </> = 0. For small wave 
vectors the excitons are superradiant with Tex larger than 
Tat- The damping rate decreases with increasing k to 
become metastable, and vanishes for k > kc- 

In figures (5) we plot T^x/^at vs. (j) iov 9 = ■n/A, 
where the transition dipole has equal in-plane and nor- 
mal components. Figure (5a) is for small wave vectors 
Eq =0.1 eV. Here the excitons are superradiant with 
maximum damping rate at = tt. In figure (5b) we plot 
for Eq = 0.9 eV. 




FIG. 4: The scaled damping rate Vex/Tat vs. So(fc) = hck, 
for e = 7r/2 and </- = 0. 
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FIG. 5: The scaled damping rate Tex/Tat vs. 0, for 9 = 7r/4 
with (a) Eq = Q.I eV, and (b) Eo = 0.9 eV. 



In figures (6) we plot Tex/Tat vs. (j) ioi 9 = 7r/2, where 
the transition dipole has only in-plane component. Fig- 
ure (6a) is for small wave vectors Eq = 0.1 eV and figure 
(6b) for Eo = 0.9 eV. The excitons are superradiant. 
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FIG. 6: The scaled damping rate Tex/Tat vs. 0, for 9 = ■w/2 
with (a) Eo = 0.1 eV, and (b) Eq = 0.9 eV. 

In figures (7) we plot T^x/Tat vs. 9. In figure (7a) we 
use (j) = Tr/2 and Eq = 0.1 eV, and in figure (7b) we use 
= ■jt/2 and Eq = 0.9 eV, where here the excitons are 
superradiant with maximum damping rate at ^ = 7r/2. 
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FIG. 7: The scaled damping rate Tex/Tat vs. 9, for 4> = 7r/2 
with (a) Eo = 0.1 eV, and (b) Eo = 0.9 eV. 

Now we calculate the emission pattern for a long wave 
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length exciton of wave vector k into free space. Using 
the Hamiltonian H = H^.^ + H^ad + Hj, we derive the 
radiation field operator equation of motion i—^^ — 
^ril) bqx(t) ~ igqx which has the formal solution 



6qA(0 = ^'qA(0)e-^"'-(«)*-gqA / dt' B^{t')e 

JO 



(6) 

In using the source term, the positive part of the electric 
field for a fixed k, at position r and time t, is 



Je(k)t 



9. A 



2€oV 



dt' Bi^(t') e 



i[iJ,{k)-LJ^(q)](t~t') 



(7) 



with r — p + zz, where we transferred the exciton opera- 
tor into rotating frame by using Bi^{t) = Bi^{t) e~*"<= 
with Eex{i^) — ?icje(k). The sum over Qz casts as before 
into an integral, and we use also the summation over the 
photon polarization J2\ GqAGqA — ^ — 

At this point we take the approximation of long wave 
length excitons, that is small wave vectors with fca <C 1, 
where k <^ Qz to obtain 



Jlii-p+q^z- 



dqzQz fJ' 



q(q - 



Aneoa'^VN Jo ' ' V Iz 

' dt' Bk(i') e*["'=('''-'=«^](*-*'). (8) 



We use = {fip^z), and q = (k, gx), and sub- 
stitute uj — cQz and duj — c dqz- Then we apply 
the Weisskopf-Wigner Approximation 12], in replac- 
ing Lu under the integrand by We(k) and to take them 
out of the integral, and to extend the lowest limit of 
the integral over to to —oo. After applying the result 
/+^ doj e*K(k)-c.](t-t'-f ) ^ 27r^ f), in terms 

of the angle (j) and in the laboratory frame, we have 



Eoik) 



2eohca^VN 
- Eo{k) 



EeA^) 



/111 COS (/) + fl. 



, -/111 COS0 I ,(9) 



where we defined the unit vector k = k/fc. 

The expectation value for the exciton operators is 

{Blit - z/c)Sk(< - z/c)) = {Blm^m e-^^'-'-^/^K 

(10) 

where Fk is the previous exciton damping rate. 

For a single exciton, where (i3^(0)i3k(0)) — 1, and in 
terms of the angle 9, we obtain the emission intensity 



E'o{k)^,' 

{2enhca^)^N 



cos 9 + sm^ 



El[k) 



El{k)-EU\,) 



■ Ei{k)-Elik) 
^ Eecc{k)Eo{k) 



+ 2 sin 9 cos 9 cos (j) 



The dependence on the distance from the lattice plane 
appears only in the damping exponent. Note that the 
result is in the limit of fca ^ 1. 



In summary our calculations predict that the radia- 
tive damping rate and emission pattern into free space of 
an exciton, for two dimensional optical lattice ultracold 
atoms in the Mott insulator phase, is strongly deviated 
from a single excited atom. The emitted light properties 
are a complicated function of the exciton wave number 
and polarization direction. For some wave numbers and 
polarizations excitons exhibit damping rate much larger 
than a single excited atom and can be considered as su- 
perradiant excitons, while for others the damping rate 
is much smaller. We found that beyond a critical wave 
number excitons become metastable with zero damping 
rate, where beyond hckc — Eex(k.c) no free space photon 
can conserve both momentum and energy with an exci- 
ton with a fixed wave vector. Metastable excitons can 
store photons for a long time and transport them over a 
long range. 



The big challenge is how to excite metastable excitons 
in such a system. As they decouple from the free space 
radiation field, they can not be excited directly by uni- 
formly shining a light on the optical lattice. However, by 
localized excitation or targeting the edge of the lattice we 
can excite at least a partial amplitude in the metastable 
domain, which should survive much longer. As an al- 
ternative one could use Raman excitation in a properly 
phase matched geometry to excite long lived excitons, 
which can be read out by a Stokes pulse. Note that some 
metastable excitons are bright states, and hence by using 
different optical elements one can excite them directly, or 
for a finite system to shine the optical lattice from the 
side. Excitons are collective electronic excitations where 
their momentum is distributed over all the lattice sites 
and not on a single atom site, therefore shining light di- 
rectly on the optical lattice in order to excite excitons 
will not destroy the Mott insulator phase. The results 
are of importance for quantum optics and quantum infor- 
mation processing, and are valid for any planar periodic 
structure of optically active materials. 
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